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Qh' Abstract 

ryK , We show that the total number of cohisions in the exchangeable coalescent process 

I driven by the beta(l, b) measure converges in distribution to a 1-stable law, as the 

pq ' initial number of particles goes to infinity. The stable limit law is also shown for the 
total branch length of the coalescent tree. These results were known previously for 

fm^ \ the instance 6 = 1, which corresponds to the Bolthausen-Sznitman coalescent. The 

n , ' approach we take is based on estimating the quality of a renewal approximation 
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to the coalescent in terms of a suitable Wasserstein distance. Application of the 
method to beta (a, 6)-coalescents with < a < 1 leads to a simplified derivation 
of the known (2 — a)-stable limit. We furthermore derive asymptotic expansions 
for the moments of the number of collisions and of the total branch length for the 
beta(l, 6)-coalescent by exploiting the method of sequential approximations. 
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1 Introduction 



H ■ Pitman [28] and Sagitov [29] introduced exchangeable coalescent processes with multi- 

" " pie collisions, also known as A-coalescents. A counting process associated with the A- 

coalescent is a Markov chain n„ = (n„(t)) with right-continuous paths, which starts 
with n particles n„(0) = n and terminates when a sole particle remains. The particles 
merge according to the rule: for each t > when the number of particles is n„(t) = m > 1, 
each k tuple of them is merging in one particle at probability rate 

Xm,k= [ x''{l-x)"'-''x-^A{dx), 2<k<m, (1) 
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where A is a given finite measure on the unit interval. The event of merging of two or more 
particles is called collision. By every collision n„ jumps to a smaller value. When A is a 
Dirac mass at the A-coalescent is the classical Kingman coalescent [23], in which every 
pair of particles is merging at the unit rate and only binary mergers are possible. Another 
eminent instance, known as the Bolthausen-Sznitman coalescent [S], appears when A is 
the Lebesgue measure on [0, 1]. 

The subclass of beta- coales cents are the processes driven by some beta measure on 
[0, 1] with density 

A(dx)/dx = --—-x''-\l-xf-\ a,b>0, (2) 

B(a, b) 

where B(-, ■) denotes Euler's beta function. This class is amenable to analysis due to the 
fact that the transition rates ([1]) can be expressed in terms of B(-, ■). For this reason and 
due to multiple connections with Levy processes and random trees, the beta coalescents 
were the subject of intensive research [21 [3], El El [9l [IHl [191 EH] • We refer to [1] for a survey 
and further references. 

In this paper we study beta-coalescents with parameter < a < 1. Specifically, 
we are interested in the total number of collisions X„ and the total branch length of 
the coalescent tree Ln- Note that X„ is equal to the total number of particles born by 
collisions, and L„ is the cumulative lifetime of all particles from the start of the process 
to its termination. The variable L„ is closely related to the number of segregating sites 
Mn, the connection being that given L„ the distribution of Mn is Poisson with mean rLn 
for some fixed mutation rate r > 0. 

Our first main new contribution is the proof of a 1-stable limit law for X„ and L„ 
as n — )■ oo. As in much of the previous work (see, for instance, [11] and [20]) we use a 
renewal approximation to n„. A novel element in this context is estimating the quality 
of approximation in terms of a Wasserstein distance. 

The second new contribution are asymptotic expansions for the moments of Xn, Ln 
and Mn for the beta (1, 6)-coalescent with arbitrary parameter 6 > 0. These expansions 
are obtained independently from the weak limiting results mentioned before. The proofs 
are based on the method of sequential approximations similar to those used in [IS] . 

The rest of the paper is organized as follows. Section [2] gives a summary of some 
results on limit laws related to the beta-coalescents. In Section [3] general properties of 
the block-counting Markov chain and basic recurrences are discussed and the main results 
are stated. Section H] recalls the definition and properties of a Wasserstein distance. In 
Section O we provide proofs of the main results. Some auxiliary lemmas are collected in 
the appendix. 

2 A summary of limit laws for beta-coalescents 

The tables in this section summarize the limit laws for X„, L„ and the absorption time 
of the coalescent r„ := min{t : n„(t) = 1}. The distributions which appear in the tables 
will be denoted as follows 

(i) A/", standard normal. 



(ii) Sa with 1 < a < 2, (spectrally negative) a-stable distribution with characteristic 
function 

2; H- exp{|2:|"(cos(7ra/2) + isin(7ra/2)sgn(2:))}, zeR, (3) 

(iii) Si (spectrally negative) 1-stable distribution with characteristic function 

2; I-)- exp { — |z|(7r/2 — ilog |z| sgn(2;))}, z G M, (4) 

(iv) £^{a, h) with a,b,'y > 0, distribution of the exponential functional J^ exp{—'ySa^b{t))dt, 
where {Sa^bit))^^^ is a drift-free subordinator with the Laplace exponent 



'^a,b{z) 



[ {l-il-xy)x''-%l-xf-^dx, z>0, 
Jo 



(v) Q, Gumbel with distribution function x i-)- exp ( — e ^'),a; G M, 
(vi) p, convolution of infinitely many exponential laws with rates i{i — l)/2, i > 2. 



Table 1: Limit distributions of (X„ — an)/bn for beta (a, 6)-coalescents. 
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Notation and comments: ri = C(2,6), r2 = 2^(3, 6), where C(')') is the Hurwitz 
zeta function; mi = ^(a - 2 + 6) - \E'(6), 7112 = \&'(6) - ^'(a -2 + 6), where \E'(-) is the 
logarithmic derivative of the gamma function. 

For the Bolthausen-Sznitman coalescent the limit law of X„ was first obtained in [9] 
using singularity analysis of generating functions. A probabilistic proof of this result 
appeared in p9], where a coupling with a random walk with barrier was exploited, and 
the technique was further extended in [20J to study collisions in the beta (a, l)-coalescents 
with a G (0, 2). The aforementioned limit laws for a > 1 are specializations of results for 
more general A-coalescents with dust component, i.e., those driven by the measures A such 
that Jq a;~^A(da:;) < 00 [IJl [131 HH [IS]. For Kingman's coalescent we have X„ = n — 1 for 
ah n eN. 

In the next two tables the value a = corresponds to Kingman's coalescent. 



Table 2: Limit distributions of [vn — a„)/6n for beta (a, 6)-coalescents. 
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Notation and comments: The constants m and s^ are 
a + b- 1 



Yl - (a + 6 - 2) (^(a + 6 - 1) - ^(6)) j , 



(a-l)(2-a) 

- (a + 6 - 2) ((^(a + 6 - 1) - ^(6))^ + ^\h) - ^'(a + 6 - 1)) J , 

Ci = h{h + 1)C(2, &), C2 = 26(6 + 1)C(3, h). The constants mi and m2 are the same as in 
Table 1, and for a > 2 

_ (a-l + 6)(a-2 + 6) 

^ ~ (a-l)(a-2) ■ 

In the case a G (0, 1), 6 > the beta (a, 6)-coalescent has the property of coming down 
from infinity [30], which implies that r„ weakly converges without any normalization to 
some limiting law, which is not known explicitly. The result for a > 1 is a special case of 
Theorem 4.3 in [I2j. The case a = 1 and 6 7^ 1 is open; in this case the coalescent does 
not come down from infinity. 

Table 3: Limit distributions of (L„ — an)/bn for beta (a, 6)-coalescents. 
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Notation and comments: The constants are a = 2 — a, P = l + a — a"^, c\ — ^_^ , 

r(a+l)(a-l)l+°~^ 



C2 



cos(7ra/2)ra ^(2-a)' 



In [21] the weak convergence of properly normalized L„ was proved for A-coalescents 
with dust component. In particular, that result covered the beta (a, 6)-coalescents with 
a > 1. Although some partial results for a G (0, 1) and 6 > were obtained in [7], this 
case with h ^2 — a remains open. 

3 Main results 

For the general A-coalescent, the Markov chain n„ is a pure-death process which jumps 
from state m to m — k + 1 aX rate (™) Am.fc, where Am,fc, 2 < A; < -m, is given by ([T]). The 
total transition rate from state m > 2 is 

^-•=E(T)^-''^= / {^-mx{l-xr-^~{l-xr)x-^k{dx). (5) 

The first decrement /„ of n„ has distribution 

The strong Markov property of the coalescent entails the distributional recurrences 

Xi = 0, X„^l + X_,„, nGN\{l}; (6) 

ri = 0, r„^T„ + <_,„, nGN\{l}; (7) 

Li=0, L„^nT„ + L'„_,„, n G N\{1}, (8) 

where T„ denotes the time of the first collision, hence T„ has the exponential law with 
parameter A„; X^ (respectively, r^, L'^) is independent of /„ (respectively, (T„,/„)) and is 
distributed like X^ (respectively, r^, L^.), for each A; G N. 
Letting A be defined by ([2]) with a E (0, 1] denote by 

Pnl-=Wn = n-k}, fc = l,...,n-l. (9) 

Using the leading terms of asymptotic relations fl27j) . fl28|) and fl29|) we infer 
hmxj^'^) (2-«)r(/c + a-l) (a) 

hence 

/„ 4 ^, n ^ oo, (10) 

where ^ is a random variable with distribution {p^ )k&i- 

Consider a zero-delayed random walk (5'„)^ defined by 



where (^j) are independent copies of ^ with distribution {p^^')keN, and let {Nn) -^ be 
the associated first-passage time sequence defined by 

Nr, = mi{k >0:Sk>n}, nE Nq. 

It is plain that 

iVo = o, iv„^i + iv;_^^„ = i + iv;_^ii^<„i, nen, (ii) 

where N'^, is independent of ^ and distributed like N^-, for each A; G N. Comparing ([6]) and 
(TTTj) one can expect that if A^„ (properly centered and normalized) converges weakly to 
some proper and non degenerate probability law then the same is true for X^ (with the 
same centering and normalization). This is what we mean by a renewal approximation 
mentioned in the Introduction. This idea was exploited in [H] (for a G (0, 1), 6 > 0) and 
in [20] (for a G (0, 1], 6 = 1) to derive the limit distribution of X„ from that of Nn- We 
shall use a method based on probability metrics to show the stable limits for a G (0, 1] 
and 6 > 0. 

Theorem 3.1. As n ^ oo the number of collisions X„ in the beta (a, b)-coalescent satis- 
fies 



(i) /or < a < 1 and b > 



(ii) for a = 1 and b > 0, 



Xn- {I- a)n d ^ 

(1 - a)nV{2-a) ^ '^2-a, 



Xn — log 72 — log log n -^ Oi. 

n 

As a consequence of our main theorem we also obtain a weak limit for the total branch 
length Ln and the number of segregating sites M„ (see [24j) of the beta(l, 6)-coalescent. 

Corollary 3.2. For the total branch length Ln in the beta. {l,b)-coales cent we have as 
ra — 7- oo 



b log^ n 



n 



Ln — log n — log log n — ;■ 5i . 



Corollary 3.3. For the number of segregating sites Mn in the beta {l,b)-coalescent we 
have as n —7- oo 

Mn — iogn — log log n — )■ Oi, 

rn 

where r > is the rate of the homogeneous Poisson process on branches of the coalescent 

tree. 

We now turn to the moments of X„, L„ and M„. An analysis of these moments 
provides further insight into the structure of these functionals. Our next result concerns 
the asymptotics of the moments of the number of collisions Xn in the beta (1, 6)-coalescent. 

6 



Theorem 3.4. Fix b G (0, oo) and j E Nq. The jth moment of the number of collisions 
Xn in the beta {l,b)-coales cent has the asymptotic expansion 



w I . rrij 



where the sequence {mj)ji^'^^ is recursively defined via rriQ := and rrij := mj_i + Hj/j for 
j e N, with Kj := ij + b- l)^(j + b)+j - {b- 1)^(6), j G Nq. 

For some more information on the coefficients rrij, j G N, we refer the reader to Eq. ( 123|) 
in the proof of the fohowing Corollary 13. 5[ which provides asymptotic expansions for the 
central moments of X„ in the beta (1, 6)-coalescent. 

Corollary 3.5. Fix b G (0, oo) and j E N\ {1}. The jth central moment of the number 
of collisions Xn in the beta. {l,b)-coalescent has the asymptotic expansion 

( — iV n^ f n? \ 

E(X„-EX„)^ = L_^B(6,j-l)--5:^ + -— ^^ , n ^ oo. (13) 

3 log"* n \\og^^ n) 

In particular, Var(X„) = (^b^~^r? j \o^ n + 0(n^/log^n) as n -^ oo. 

Remark 3.6. For & = 1, Eq. flT3|) reduces to the asymptotic expansion (see Panholzer [27t 
p. 277 or Theorem 2.1. with a = 0]) 

E(X„-EX„)^ = . . \.^ ,+1 +0^-^T^ , n^oo 
j(j -l)log-^n xXo-g n) 

of the jth central moment of the number of collisions X„ for the Bolthausen-Sznitman 
n-coalescent. 

The last result concerns the moments end central moments of the total branch length 
Ln of the beta (1, 6)-coalescent. 

Proposition 3.7. Fix b G (0, oo) and j G Nq. The jth moment of the total branch length 
Ln of the beta (1, 6)-coa/escent has the asymptotic expansion 

1 Tt^ I m ■ I 1 

^Un = - — — l + -^ + — ^11, n^oo, (14) 

fe^log^nV logn \\oin ' '' ' ^ ^ 



where the sequence (mj)jgNo is defined as in Theorem \3.4\ Moreover, for j G {2,3, . . .}, 



the jth central moment of Ln has the asymptotic expansion 



jb^ log-' n Vlog n 



E(L„-EL„)^^ = i-^B(6,j-l)-^— ^^ + -— ^2-), n^oo. (15) 



In particular, Var(L„) = {2b^) ^n'^ / log"^ n + 0(n^/ log'^n) as n ^i- oo. 

Proposition 13.71 indicates that 6L„ essentially behaves like X„, in agreement when 
comparing Theorem 3.1 (ii) with Corollary 3.2. The proof of Proposition 13 . 71 works essen- 
tially the same as the analogous proofs of Theorem 13.41 and Corollary 13.51 for X„. Instead 



of the distributional recurrence (jH]) for (X„)„gN one has to work with the distributional 
recurrence ([8]) for {Ln)n&N- Since the expansion of ET„ = 1/A„ is known (see Lemma 
I6.4p . the proofs concerning X„ are readily adapted for L„. A proof of Proposition 13.71 is 
therefore omitted. We finally mention that, for the beta (1, 6)-coalescent with mutation 
rate r > 0, expansions for the moments and central moments of the number of segregat- 
ing sites M„ can be easily obtained, since (see, for example, [HI p. 1417]) the descending 
factorial moments of M„ are related to the moments of L„ via E(M„)j = HEL^, j G Nq. 

4 Probability distances xt and dq 

For real-valued random variables X and Y and T > the Xr-distance between X and Y 

is defined by 

Xt{X,Y) = sup |Ee^*^ - Ee^*^|. (16) 

\t\<T 

By the continuity theorem for the characteristic functions convergence in distribution 
Zn — )■ Z holds if and only if lim xriZn, Z) = 0, for every T > 0. 

n— >oo 

Let Vg, q G (0, 1], be the set of probability laws on M with finite gth absolute moment. 
Recall that \x — y\'^ is a metric on M. The associated Wasserstein distance on Pg is defined 

by 

dg{X,Y)=miE\X-Y\'^, (17) 

where the infimum is taken over all couplings (X, Y) such that X = X and Y = Y. 
For ease of reference we summarize properties of dg in the following proposition. 

Proposition 4.1. Let X,Y be random variables with finite qth absolute moments. The 
Wasserstein distance dg has the following properties: 

(Dist) dq{X,Y) only depends on marginal distributions of X and Y , 

(Inf ) the infimum in (I17p is attained for some coupling, 

(Rep) the Kantorovich- Rubinstein representation holds 

dg{X,Y)= sup|E/(X)-E/(F)|, 

where 7g := {/ G C(M) : \f{x) - f{y)\ <\x - y\\ z, y G M}, 

(Hom) dg{cX, cY) = \c\id{X, Y) for ceR, 

(Reg) for X,Y,Z defined on the same probability space dg{X + Z,Y + Z) < dq{X,Y) 
provided Z &T>g is independent of (X, Y), 

(Aff) dg{X + a,Y + a)= dg{X, Y) for aeR, 

(Conv) for X, X„ G Vg convergence dg{Xn,X) — )■ 0, n — ?■ oo implies X„ — > X and 

EIXJ^^EIXK. 



Proof. We refer to [TT| ^T\ for most of these facts. To prove (Reg) choose an independent 
of Z couphng (X', Y') on which the infimum in the definition of dq is attained. Then 

X + Z = X' + Z,Y + Z = Y' + Z and the definition of dg entails 

dg{X + Z,Y + Z) < E\{X' + Z)-{Y' + Z)\'^ = E\X' -Y'l'^ = dg{X,Y). 

Property (Conv): the convergence of moments is easy; the rest is a consequence of Lemma 
gj to follow. D 

Lemma 4.2. For T > and q G (0, 1] there exists constant C = Cx^g > such that 

sup |Ee**^ - Ee**^| < Cdg{X, Y), neN. 

\t\<T 

Proof. Assume that the infimum in the definition of dg{X, Y) is attained on {X, Y). It is 
easy to check that for arbitrary q G (0, 1] 



e'^ -e'y =2 



X — y 
sm 



< 2^-HIg\x-y\'^, x,y e 



2 
where Mg := sup^^Q | smu\u^'^ < oo. Hence 

sup |Ee**^ - Ee^*^| = sup |Ee**^ - Ee**^| < sup E|e^*^ - Ee'^^\ 

\t\<T \t\<T \t\<T 

< 2^-'iMg sup |t|'?E|X -Y\'^ < 2^-'^MgT''dg{X, Y), 

\t\<T 

as wanted. D 



5 Proofs 

5.1 Proof of Theorem 13.11 

Suppose a = 1. It is enough to show that 

Tog^n 



hm XT { Xn - log n - log log n,Si) = 0, 



for every T > 0. 

Using the triangle inequality yields 



XT y ^n - log n - log log n, 5i j < 

/log^n log^^ 

Xt Xn — log n — log log n, Nn — log n — log log n 

\ n n 

,2 , 



/ log n \ 

+ Xt{ A^„ -logra-loglogra,5ij. 



The second term converges to zero by Proposition 2 in [12] on stable limit for the number 
of renewals. In view of Lemma 14.21 to prove convergence to zero of the first term it is 
sufficient to check that 

/ lo£C Ti iOS. Tl \ 

lim dq ( Xn — log n — log log n, iV„ — log n — log log n ) = 0, 

n-s>oo \ n n J 

for some q G (0, 1], which in view of the properties (Hom) and (Aff) in Proposition 14.11 
amounts to the estimate 

cig(X„,Ar„) = o(nMog-''^n), n ^ oo. (19) 

In the like way, proving Theorem 13.11 in the case a G (0, 1) reduces to showing that 

rf,(X„,iV0=oK/(2-a))^ ^^^^ ^20) 

for some q G (0, 1]. 

Using recurrences ([6]) for X„ and flTTl) for A^„ we obtain 

K := rf,(x„, iv„) = d,(x;_,^, iv;-(CA„)) < ^.(K-/„, K-(ca„)) + ^^^(x;.,^, iv;_,j 

n-\ 

< dq{K_j^, K-i^An)) + E|X„-/„ - Nn-iS =: Cn + J] P{/„ = n - A;}E|Xfc - iV,|^ 

for arbitrary pairs (^{Xk, N^)) -^^^f^^^_^ independent of /„ such that Xk = Xk, Nk = Nk. 
Passing to the infimum over alFsuch pairs leads to 



n-l 



tn<Cn + J2 ^{^n = n- kjt^. (21) 



fc=l 



We shall use (1?I1) to estimate t„. 

First we find an appropriate bound for c„. Let (/„,0 be a coupling of /„ and C, such 
that (recall (Inf) in Proposition 14. ip dg{In,^ An) = E|/„ — ^Anj^. Let {Nk)i^^ be a copy 

of (A'^fc)^ J5J independent of (/n,0- Since [in,^, i^k)) is a particular coupling we have 

Exploiting the stochastic inequality 

N,+y-N,<Ny, x,yeN 



yields 



EIAT f - N f I"? < EX''- . 



Furthermore, we obviously have X„ < n, hence 

c„ < E|/„, - ^ A n|'' = (ig(/„, ^ An). 

10 



Now we invoke the Kantorovich- Rubinstein representation ((Rep) in Proposition 14. ip for 
dq. Set J-'q^o := J-'g n {/ : /(O) = 0} and note that / G J^q^ imphes \f{x)\ < |x|'', x G M. 
We have 

c„ < dqiln, CAn) = sup E/(/„) - E/(e An) = sup |e/(4) - Ef{C A n] 

n — 1 n—1 

= sup |5^P{J„ = A;}/(fc)-^P{e = fc}/(fc)-/(n)^P{e = A;} 

/e^g,0 ^i.^]^ ^^]^ ;,>„ 

n-1 



< ^\¥{In = k}-¥{i = k} k" + n''W{i>n} 



k=l 

For appropriate g G (0, 1] (to be specified below) such that a + q > 1 use Lemma 
16.31 in the Appendix along with the relation F{C, > n} = 0{n°'^^) to obtain the estimate 
Cn = 0(n''+"~^). With this bound for c„ a 0-estimate for t„ follows using Lemma [Ol 

If a G (0, 1) one can take q = 1. Then the cited lemma applies with ipn = n and 
Tn = Mn°'~^ (M large enough) and gives estimate 

dq{X^,K,) = 0{n-), 

which implies fl20l) . 

For the case a = 1 application of the same lemma with ijjn = n/{log{n + 1)) and r„ = 
Mn'^-i (M large enough) leads to t„ < Mra'?(logra)-^ Thus ([H]) holds for q G (0, 1/2). 
The proof is complete. 

5.2 Proof of Corollaries [32] and [3731 

We follow closely the proofs of Theorem 5.2 and Corollary 6.2 in ^. In view of 

!^i!^L„ - log „ - log log n = !!^A-„ - log „ - log log „ + l^^ (te„ - A',) , 

n n n \ / 

it is enough to show that °^ " ( 6L„ — X„ j — )■ in L2. 

Let Tj's be independent exponential variables with rates Aj, j > 2. Assuming the T/s 
independent of the sequence of states visited by n„ we may identify Tj with the time n„ 
spends in the state j provided this state is visited. Given the sequence of visited states is 
n = -io > "^1 > ■ ■ ■ > ik-i > ik = ^ the total branch length L„ is distributed like J2r=o ^rTtr 
for n G N\{1}. 

For k G {1, . . . , n} and i = (29, ■ ■ ■ ,ik) with n = io > "^i > ■ ■ • > "^fe-i > "^fe = 1 define 
the events A^^i := {X^ = k, (n„(to), • • • , n„(tfc)) = i}, where to = and ti < ^2 < • • • are 
the collision epochs. We have 

fc-i 2 

EibL^-X^f = 5^P{Afc,i}E(5^(6vT,, -1)) 

k,i r=0 

k,i r=0 r,s=0,ry^s 
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Furthermore, A„ = 6ra + 0{logn) as ra — )■ oo for a = 1 and b > (see (1291) ) which imphes 

\E{bkTk - 1)1 = 0{k-Hogk) and E{bkTk - if = l + 0{k-Hogk). 
Therefore, 

n n r, 

E(6L„-X„)2 < J]P{Afc,i}(j]E(6rT,,-l)2+(^|E(6rT,-l)|) ) 

fc,i r=2 r=2 

= ^P{Afc,i}(^n + 0(log^n)) =n + 0(log^n), 

A;,i 

and the convergence in L2 follows. 

Corollary 13.31 follows from the fact that given L„ the distribution of M„ is Poisson 
with mean rL„. See Corollary 6.2 in ^ for details. 

5.3 Proofs of Theorem 13.41 and Corollary 13.51 

Let us verify (TT2l) by induction on j G N. From ([6]) it follows that ai := EXi = and 
a„ := EX„ = 1 + X]m=2P".'»^™' ^ £ N \ {1}. In the following we apply the method 
of sequential approximations to the sequence {an)nm- The sequence (6n)neN, defined via 
61 := and bn := an — n/ logn for n G N \ {1}, satisfies the recursion 

n— 1 ^ \ n— 1 

bn = an- = 1 + y] Pnli { 7^^ + ^m ) " r^^ = ?« + V P^^L^m, 

logn ;^2 Vlog"^ / log^ ^2 

n G N\ {!}, where Qn := I — n/\ogn + X]m=2P":'""^/l'^g"^' n E 'N\ {!}. By Corollary 
16.71 (applied with a := 1 and p := 1), 

logn \logra logn \\og nj J logn ylog ""-/ 

where mi := Cb_i,i = 2 + \1'(6). The sequence (c„)„gN, defined via Ci := and Cn := 
bn — rriin/ log nfornGN\{l}, therefore satisfies the recursion 

n— 1 ^ \ n—1 

Cn = bn-mi-\- = qn + J2Pn!m(^l7^^ + Cm]-mi-^ = % + J^ Pn^rnCm, 

log n ^^ ' \ log m J log n ^^ 

n e NX {!}, where q'^ := Qn - min/log^n + mi Z]m='2 P",™"^/ log^ "^^ ra G N \ {!}. By 
Corollary 16.71 (applied with a := 1 and p := 2), 



n 



Qn = Qn-mi- 2— +"^l(^ 2 ^i ^ '^ \ \ 2—)) = ^i- ^ (' 

log n Vlog IT- logn \log nJ J ylog n/ 

since g„ = mi /logn + 0(l/log^n). By Lemma 16721 (applied with a := 1 and p := 3), it 
follows that Cn = 0{n/ log n). Thus, (IT2l) holds for j = 1. Assume now that j > 2. From 
EX£ = E(X„ - 1)^' = Y^rl 0(-l)J-*EX; + EX^ it follows that 

j-l ^ .^ n-l 

an, := EX^ = 5^(]](-l)^-^-EX; + EXi;^ = g„_^. + ^pa)^a^,^., 

2=0 ^ ^ m=2 
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neN\ {1}, where qnj := Eto (i)(-l)^"^"'E-^n' neN\ {1}. Since, by induction, for 
all i < j, 



log^raV" ' log"- ' "'Vlog^^, 
it follows that (the sunimand for i = j — 1 asymptotically dominates the others) 

log-' n V log""- Vlog^ 

Now apply the method of sequential approximations to the sequence (a„j)neN- The se- 
quence (&n,j)neN5 defined via bij := and 6„j := a„j — n^ / log' n for n G {2,3,...}, 
satisfies the recursion 



n-l 

E 

m=2 



Vi = 9n,j + Y^ Pninbm,j, n G {2, 3, . . .}, 



where g^^^- := g„j- - n-^'/ log-^ n + Y,m=2 Pnlnm^ / log-' m, n e {2, 3, . . .}. By Corollary [631 
(applied with a := j and p := j). 



log-^ ra log-^ n \log-^ nj 



n J log-' n 

H ^ j ^-j h K j ■■ h O r— — 

log-' n log-' n log-' n \log-' n 

n^~^ ( n^^^ 



log-' n ylog-'''"^ 



n 



where Kj := Cbjj and rrij := wij_i + K^/j. The sequence {cnj)neN, defined via Cij := 
and Cnj := bnj — rUjU^ / log-'^^ ra for ra G {2, 3, . . .}, therefore satisfies the recursion 



Cn,j — Qnj ' / J Pn,m'^rn,jy fl t \^; -J; • • -J; 



n-l 

E. 

m=2 



where q'^j := g^ ,,• - mj-n-'/log-'^^ n + m^ E^J2Pn,m'm-'/log-'"^^m, n G {2,3,...}. By 
Corollary 16.71 (applied with a := j and p '■= j + 1), 

<ln,i = Jmj-—- + O -—rr- - ^j 



log-' n \\og^^^ nJ log-'"''^ 



n 



+mA-^T^-3-^ + 0\-^-^\ = O 



\ loe;-''^ n log-' n \ log-''^ nil \ log-''*' 



n 



.i+2^N 



By Lemma [6l2] (applied with a := j and J9 := j + 2), it follows that c„j = 0[p? j log-' 
which shows that (TT2|) holds for j. The induction is complete which finishes the proof of 
Theorem 13.41 
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We now turn to the proof of Corollary 13.51 Let us first verify that the sequence 
(?7ij)jgNp, recursively defined in Theorem 13 ■4[ satisfies the inversion formula 

J2(l){-iy-'m, = t]^B{b,j-l), JEN\{1}. (22) 

i=o ^ ^ -^ 

Using the formula \E'(a; + l) = \E'(x) + l/a;, x G (0, oo), it is readily checked that Kj^i — Kj = 
2 + ^(6 + j), J e No- For all j E Nq it follows that Kj = Etol'^i+i - f^i) = Eto(2 + 
^ib + t)) = 2j + J2lZ^ ^{b + t) and 

^j = Y.{m-m,_,) = E f = E (2+7 E ^(^+^)) = 2j+E ^(^+^) E T (23) 

By O, for J 6 {2,3,...}, 

E(->-i)'-'"< = tC)(-ir*(2'+E*(''+*) E }) 

j=0 ^ ^ i=l ^ / V fc=0 Z=A;+1 ^ 

i=l ^ ^ A:=0 i=fe+l 

£*("+«=) E !E (■)(-!)'" = £*(''+*) E KrI)<-^>'"' 

fc=0 i=fc+l i=l ^ ^ fc=0 i=fc+l ^ ^ 

ig*(6+*) t 0')(-ir' = 7 £*(''+ *)(';')'-'' 



fc=0 i=fc+l j=« ^ ^ fc=0 i=fc+l 
J-1 J' /-A 1 ^-1 

•^ A;=0 l=k+l ^ ^ -^ fc=0 

Plugging in p^^) = (;^I^) + ("'^^) and reordering with respect to p^^) leads to 



^(^V-ir^m, = ^X^(-ir^-^K^)(v&(6 + fc + l)-v&(6 + A:)) 



i=0 ^ ^ -^ fc=0 

J-2 



? 

•^ fe=0 



where the last equality holds, since Y12=o(~^)^ (T) / ^^ + k) = B{b,n + 1) for all n e Nq, 
which is for example readily verified by induction on n G No- Thus, (!22l) is established. 

Thanks to Theorem 13.41 and the inversion formula fl22l) the proof of Corollary l3.5l is now 
straightforward. Basically the same argument has for example been used by Panholzer 
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[271 P- 277]. Plugging in the expansion flT^ for the ordinary moments shows that 
■^ Vv log'nV log'^ \\og^nJ J \\ogn\ logn Vlog^ri 



log-'np^Vv V logn \\ognJ J \ logn \\og n 



E Mr- iH-M^^^i^+o 



log-'n'^Vv V logn \log n 



? C) <-^''"" + ii ? (0 <-^''"'«^' - ''"■ + "'' + « ( 



log^ n ir'r, Vv W^^ nj-^\ij ' Vlog^^^ 



i=0 ^ ^ *= j=0 



n 



log-'+^n j ' Vlog-^"^^ 



n 



since, for J > 2, ELo Qi-^Y" = 0, Eio 0(-l)^-^(j-O = 0, and ELo (9(-l)^-^^. = 
{—iy/jB{b,j — 1) by (122!) . The proof of Corollary 13.51 is complete. 

6 Appendix 

For each n G N let {pn,k)o<k<n be an arbitrary probability distribution with p^.n < 1- 
Define a sequence (a„)„gN as a (unique) solution to the recursion 



n 
k=0 



n + ^Pn.feOfc, nen, (24) 



with given r„ > and given initial value ao = a > 0. The following result is Lemma 6.1 
from [12j. 

Lemma 6.1. Suppose there exists a sequence {ipn)n£N such that 

(CI) liminf„^oo ^n ELo(l ~ k/n)pn,k > 0, 

(C2) the sequence {rkipk/k)km is non-increasing. 

Then an, defined by ( l24l) . satisfies 

«- = 0(E^)' ^^^- (25) 

k=l 

Lemma 6.2. Let {an)n&N be a sequence of real numbers satisfying the recursion ai = 

and an = Qn + Em=2 P".™'^'" ^ n e N \ {1}, for some given sequence (gn)neN\{i}- U 
Qn = O^n""'^ / log^~^ n) for some given constants a G (0, oo) and p G [0, oo), then an = 
0{n"/logPn). 
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Proof. Fix some S such that < 6 < a. Set a'^ := \an\/n^ and q'^ := |g„|/ra'^. Then 
Qn ^ Mn"^^^^ / log^~ n =: r„ for some M > and all n > 2. Further, 



n—l I I n— 1 



"n — 'in ^ / ^ yn,m. x — ^n ^ / ^ Pn, 



n" " '■ — '" ■"'"" wP 



in=2 m=2 

n—l n—l 



ni=2 m=2 



o(^) a' 

■r n,mm' 



Set V^n := n/logn, then both conditions (CI) and (C2) of Lemma [6.11 are fulfilled. Hence 
«n = 0(^^=2 ^°"^"V log'' ^) = O in''-^ / \og^ n) and |a„| = n^a'^ = OK/log^n). D 

Lemma 6.3. For the first decrement In of the Markov chain (n„) associated with the 
beta, {a, b)-coales cent (a G (0,1] and b > 0) and a random variable ^ with distribution 

n-l 

Y^ A;''|P{/n = k}-¥{i = k}\= 0(n"+'?-2), (26) 

fc=i 

whenever < g < 1 and q + a> 1. 

Proof. For the beta (a, 6)-coalescents formula ([1]) reads 



^ Bia + k-l.n-k + b-l] 



An,fc+1= / x''-\l~xY-''-^k{dx) 

/o 



Using the known estimate for the gamma function (see formula (6.L47) in [T]) 

< Mc X-'^-\ n > 2, c, d > -2, 



r(n + c) ^,_, 



r(n + d) 



we obtain 



n \ _ / n \B{a + k~l,n-k + b-l) 

k + ir''^' - \k + l) BM^ ^ ^ 

V{n + l)r(a + A; - l)r(n - /t + 6 - 1) 



r(A; + 2)r(n - k)T{n + a + 6 - 2)B(a, 6) 
r(a + A;- 1) 



(A; + l)!B(a,6) 



(ri^-a-b _^ Q(^2-a-6)j /(^ _ ^)b-l + o('(n - fc)^-^ 



uniformly for 1 <k <n — 1 and n > 2. 

Using (|5]) with A given by ([2]) we infer (see also Corollary 2 in |14] ) 



, _ ^M n'-^ ^ Oin^-^) = ^-^ 

" (2-a)B(a,6) ^^ ' (2-a)B(a,6)' 



A„ - -^^^^^-.n^- -H 0(„-") ^ y^^^^lbr"^-" (l + 0("-')) . (2B) 



when a G (0, 1) and 6 > 0, and 

A„ = 6n + 0(logn), (29) 
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(") _ (2 - a)r(a + A: - 1) ^,_, /^^ _ ^^,_, ^ ^ / ^^ _ ^^.-2 ) w ^ ^ q(^-1) 



when a = 1 and 6 > 0. Hence for 0<a<l,6>0, 'n,>2 and k = 1, . . . ,n — 1 

^"'"-'^ ~ r(a)(A; + l)! 

= pi'^) ((1 - k/nf~' + 0(n-\l - A;/n)^-2)) (l + O(n-i)) 

= p(") ((1 - k/nf-' + 0(n-\l - A;/n)^-2)) 

= pi^\l - k/nf-' + 0(p'j^^n-\l - k/n)'-' 

Analogously for a = 1 

Pn!n-k = p'k\{l-k/n)'-' + o(n-\l-k/nf-^^^(^l + 0{n-'\ogni^ 

= p^^^ ({I - k/nf-^ ^0(n-\l- k/nf~^) + O (n'^ log n(l - k/n)^-^) ) 

= p^^\l - k/nf-^ + 0Uk'^n-\l - k/nf-^) + oUk'^n-Hognil - k/nf-'Y 

Substituting these expansions into the left-hand side of (l26ll gives 



n-l 



n-1 



J2k''mln = k}-F{^ = k} < J^Pini^ 



k\b'l 



fc=l 



fc=l 



n 



n-l 



- ^Eri"'Ki-!r 



Cl 

n 
fc=i 

=: 5*1(0, ra) + 5*2(0, ra) 



nJ 



for < a < 1, and 

"~-^ 1 ""^^ U ft— 1 

J^k'kUn = k}- P{^ = k} < 5i(l,n) + 52(1, n) + ^^^^^p«A;^[l - -) 



fc=i 



n 



k=l 



--: 5i(l,n) + 52(1,72) + 53(1, n). 



for a = 1. Here and hereafter ci,C2,... denote some positive constants whose values 
are of no importance. Our aim is to show that Si{a,n) = 0{n'^^°'~'^) for i = 1,2 and 
53(1, n) = 0{n'^~^). By virtue of p^ < c^k"'~^, for all A; G N, we infer 



5*1 (a, n) < 


C3)_^P+'^-^ 
fe=l 


( k\b'l 

1-- -1 
V n/ 










= 


[n/2] 

C3 )^ A:"+^-3 


/ k\b-i 

1 - 1 

V nJ 


+ C3 )J r+''-3 

fc=[n/2]+l 


1 - 1 

V nJ 


< 


[n/2] 

n ^-^ 

k=\ 


\n 

k 


= K2] + 1 


(^ 


77,/ 


)■ 
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where the inequahty |(1 — x^ — 1| < c^x, x G [0, 1/2] has been utihzed. The expression 
in the parentheses converges to J-^^,^x"'^'^~^\{l — x)^~^ — l|da; < oo. Hence Si{a,n) = 

Similarly 



S2{a,n) 



n ^-^ V nJ 

k=\ 



n ^-^ V n) n ^-^ V ru 

k=\ A;=[n/2] + l 

1-- 

n ^ — ' V n/ ^ — ' \ n/ 

fc=l fc=[n/2] + l 



fc=l fc=[n/2] + l 



n 

k=l k=[n/2]+l 

since the first term is 0{n'^) and the second term is 0(n"^''~^) by the same reasoning as 
for 5*1(0, n). 
Finally 

"^1 7 ;. 1 






k=l 

n—l T 1 1 I n— 1 



< 



C7logn^,„_2 



n 



k\b-i 



yk^-' 1--] -1 +^:^^^^yA; 



Cy log n 



.9-2 



A:=l fc=l 

= 0(n^~^ log n) + 0(n~"'^ log n), 

in view of the estimate for Si{a,n). Therefore Ss{l,n) = 0{n'^~^) and the proof is 
complete. D 

We provide a basic lemma concerning the total rates of the beta (1, 6)-coalescent. 

Lemma 6.4. The total rates A„, n eN, of the beta {l,b)-coales cent are explicitly given 
by 

n—l J 

A„ = by- = b(n-l]-b(b-l)(^(n + b-l)-^(b)), neN. (30) 

^^-^ b + k — 1 

fc=i 

Moreover, the total rates have the asymptotic expansion 

X^ = bn-b{b-l)\ogn-b + b{b-l)'^{b) + 0{n-^), n ^ 00, (31) 
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and the inverse of the total rate A„ has the asymptotic expansion 



X„ bn\ n n \ n 



^ = T^(l + ('-l)^+ ^' . ' ^' +0(^^)). n^^- (32) 



Proof. Eq. (15U]) is well known (see, for example, [T71 Appendix, Eq. (19)]. The expansion 
(13 ip follows directly from (15U]) . since '^{n + h — 1) = \ogn + 0{n~^) as n — )• oo. The last 
assertion f l32|) follows from 



6n logn l-(&-l)^(&) 

— - 1 - 6 - 1 

An n n 

hn^ - A„(n + (6 - 1) logn + 1 - (6 - 1)^(&)) 

nA,i 

O(log^n) flog^n 

nXn \ n? 



(33) 



where the very last equality holds, since A„ ~ bn, and the equality before follows by 
plugging in (l3Ti) for the term A„ occurring in the numerator of the fraction in (l33ll and 
multiplying everything out. D 

The next Lemma 16.51 provides an asymptotic expansion as n — )■ oo for sums of the 
form 

n-l a 

p , aeR,pe[0,oo). 

_ [n — m)[n — m + I) log m 

For parameters a > we will need an even sharper version (see Lemma [6l6] below) . but we 
start with this simpler version, which holds for arbitrary a; G M. Given the overlap with 
the proof of the following Lemma 16.61 and given the fact that the proof is considerably 
simpler than that of Lemma 16. 6[ the proof of Lemma 16.51 is omitted. 



Lemma 6.5. For a €iM. and p G [0, oo), as n —^ oo, 



n-l 



7 77 TTl 5 = \ 5~ 1 ~ '^ ^^ 

In — m ] r) — m -I- \ ] incr^ m incr^ r) \ 11 



m=2 



n — m){n — m + l)\og^ m log^n\ n \n 



log n r\f ^ 



The following Lemma 16^ is a sharper version of Lemma 16751 with the cost that it holds 
only for a > 0. It will turn out (see the following Corollary 16 . 71 and the proof of Theorem 
13. 4p that the expansion in Lemma 16.61 is fundamental for the analysis of the moments of 
the number of collisions of the beta (1, 6)-coalescent. 

Lemma 6.6. For a G (0, cxd) and p G [0, oo), 

7 77 ,,, „ = :-^^ l-a^^ + ^^ — - + 



m=2 
as n —;■ OO. 



{n — m){n — m + l)\og^ m log^'ny n n ynlogn 
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Proof. Note first that 

Y^ m°' /I — log(r7i/?2) 

■^ (n — m){n — ?Ti + 1) ylog^n log^^"^ n 

1 ^^ m p -^-^ 

o-P n ^ i T, — m\i T, — m ^ ^\ ]^rrP+^ n ^ 



p ^-^ m°'{—\og{m/n)) 



loe^ n ^-^ (n — m) (n — m + 1) loe:^''"^ n ^-^ (n — m) (n — m + 1) 



1 af, ahgn a^(a) 
-n [ 1 h 



log^n \ n n 

I 1 ^^ L^ — - + 



\nlognJ J log^^ n 



log^ n\ n n \n log n 

Thus, it suffices to verify that 

n-l 

_ {n — m){n — m + 1) \log^ m log^ n log^"*"^ n J Vlog^^"*^ ^ 

The function /„p : (l,ri] — )• M, defined via 

1 1 — log(x/n) 



v^ 1J^ [^ ^_ -log(m/n) \ ^ ^/ ^° ^ \ . . 



^"^^"'^ ■ log^x log^n ^"log^+^n 



^"p^^) = :; w^Tr--r-3;Tr-:) < 



has derivative 

p 

x\\og^^^n log^"^^a;, 

and satisfies fnpiji) = 0. Thus, fnp > 0. In order to verify (^^ we use a decomposition 
method. We spht up the sum on the left hand side in (jM]) into two parts Ylm=2 ■ ■ ■ ^"^^ 
X]m=a +!•••' ^^^ handle these two parts separately. We work with the sequence (a„)„gN 
defined via ai := 1 and a„ := [n/ log''^^ n\ for n >2. For the first part we obtain 

Em „ , , v-^ m 1 

7 ry -T-Jnv{m) < J] 
„m — mm — m + 1 -^^ 



(n — m)(ri — r/i + 1) -^-^ (n — ?7i) (n — ?7i + 1 ) log^ ?7i 



m=2 ^ ' ^ ' m=2 

~ log^ 2 ^ \n — m n — m+1 ) ^ \ log*''*'^ n J ' 

m=2 ^ / \ o / 

since 

l' 1 = = ^«"-l 

\ n — 777, n — m. -\- 1 I n — a„ 77 — 1 fr;, - 

m=2 



E^ 



n — m n — m+lj n — an n—1 (n — an){n — 1) nlog^^^n 
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Moreover, for the second part we have 



n— 1 , 

m 



Yl 7 V —r^fupim) 

m=a„+l ^ ' ^ ' 

n— 1 



Em" /I 1 — log(m/n) 

fn — m\ (n — m -I- 1 "l V locr^ m Incr^ n 



m 
n—i 

E 



-^ (n — m) (n — m + 1) Vlog^ m log^ n log^^^ n 

171°' log^ n — log'' m — p{— log{m/n)) log'' m/ log n 



— ^ (n — ■m)(n — m + 1) \og^ n\og^ m 



n— 1 



< 



n" Y^ log^n — log^m + plog(m/n) log^'m/logn 



log^ n log^ a„ ^f-^ (n — m) {n — m + 1) 

n" Y^ log^ n — log*' m + p log(m/n) log*' m/ log n 



^^^''^m^+i (n-m)(n-m + l) 

since loga„ ~ logn as n — ?■ oo. Thus, it remains to verify that 

5^-^ log^'n — log^'m + plog(m/n) log*'m/logn Aog^^^ n\ , , 

^^ (n — m)(n — m + l) \ n J 

m=an + l 

Let us distinguish the two cases p > 1 and p < 1. Suppose first that p > 1. Then the 
map a; I— 7- x** is convex on [0, oo). Thus, y^ — x^ < p{y — x)y^~^ for all x,y & [0, oo) with 
X < y. It follows that y^ — x^ + p{x — y)x^/y < p^{y — x^y^^'^ for all x,y E [0, oo) with 
X < y. Applying this inequality with < a; := logm < y := logn yields 

■^^ log^'n — log^'m + plog(m/n) log*'m/logn 

m=an + i 

^-^ log** n — log** m + p log{m/n) log** m/ log n 
~ ^-^ (n — m)(n — m + l) 

m=l ^ ' ^ ' 

n—1 9 /, , NO , T>— 2 n—1 1 2/ / \ 

>r^r(logn-logm)^log*' n 2^ p^2 sr^- log (m/n) 

< > = p log^ n y -. 

^^-^ {n — m){n — m+l) ^-^ in — m){n — m + 1) 

m=l ^ ' ^ ' m=l ^ ' ^ ' 

Note that 

■^-^ log [m/n) 1 Y^ log (m/n) 



Y^ log {m/n) 1 Y^ 

^^ in — m)(n — m + 1) n ^-^ fl 



m- 



'' (n — ni){n — m + 1) n ^-^ (1 — m/n){\ — [m — l)/n) 

where the particular value F (3)^(2) of last integral is obtained by choosing s := 2 in the 
chain of equalities 



log(l-M))*j^ /"^ x'e" j^ p sx 





2 du = — —^Ax = / ^ — -dx = r(s + l)C(s), s > 1, 

w^ ./n e^ - 1 V' ./n e^ - 1 
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which are based on the substitution x = — log(l — m), partial integration, and on formula 
23.2.7 in [1]. Thus, the expression on the left hand side in (135|) is even 0((log^~^n)/n). 
In particular, (135|) holds. Suppose now that p G [0, 1). Then the map x i— )■ x^ is concave 
on [0, oo). Thus, yP ~ x^ < p{y — x)x^~^ for all x, y G (0, oo) with x < y. It follows that 
yP — x^ + p{x — y)x^/y < p{y — x^x^^^/y for all x,y E (0, oo) with x < y. Applying this 
inequality with < x := logm < y := logn yields 

^^ log^n — log^m + plog(m/n) log^?7i/logn 
~ ^^ (n — m)(n — m + l) 

m,=a„ + l 



< 



< 



< 



^^ p(logn — logm)^(log^ "'^m)/logn 
^-^ (n — m)(n — m + 1) 

plog^~^an v;^ log^(m/'ri) 



E 



logn ^-^ (n — m)(n — m + 1) 

pXo^'^ansr^ log^(m/n) _ (\o^~'^ 

< ^ ( Ti — m iiT) — m -I- 1 1 



loffrz ^^-^ (n — m)(n — m + I) \ n 



m=l 



n 



since loga„ ~ logn and the last sum is 0(l/n) as shown above. Again, f l35|) holds. D 

The following Corollary 16.71 is essentially obtained by combining the three Lemmata 
16.41 16.51 and 16.61 It provides an asymptotic expansion for the sum Y^^=2Pn!mTn'^ / \o^ m. 
This expansion is a key tool for the proof of Theorem 13.41 

Corollary 6.7. Fix a G [1, oo) and p G [0, oo). For the beta (1, b)-coalescent with param- 
eter b G (0, oo), 

i:pi.;!„r^ = r^(l-"i^+^ + o(-^)). n^oo. (36) 
^-^ log m log^ n\ n n \n log nj J 

where Cb,a,p ■= (a+6-l)^(«+&-l)+p+l+(l-&)^(&) = (a+6-l)^(a+6)+p-(6-l)^(6). 

Remark 6.8. The following proof shows that Corollary 16. 71 even holds for the slightly larger 
range of parameters a,b E (0, oo) satisfying a + 6 — 1 > 0. However, we need Corollary 
16.71 only for a G [1, oo) and b G (0, oo), in which case a + b — 1 > automatically holds. 

Proof. Let gnm '■= A„P{/„ = n — m} denote the rate at which the block counting process 
moves from the state n to the state mG{l,...,n — 1}. It suffices to verify that 

n-l a 

Em 

m=2 ^ 

log^ n\ n n \n log n J J 

since ( l36l) then follows from pn,m = gnm/^n by multiplying ( l37l) with ( l32l) . Note that 

n,! 1 r(6 + m — 1) 



r{b + n — 1) {n — m){n — m + 1) {m — 1)! 
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1 < m < n. 



Since the first fraction lias expansion 



' '^-('-JV-^0(^)l (3. 



T{b + n — l) n''~2 \ \ 2 Jn \n 

it lience suffices to verify tliat 



n-l 



■^-^ 1 r(6 + r7i — 1) m°' 

^-^ [n — 'm){n — m + 1) [m — 1)\ log^m 



log n \ n 

2 



^-^)+{a + h-l)^{a + h-l)+p ( 1 



,0- , (39) 

n \n log n J J 

since ( 137|) then follows by multiplying ( 139|) with ( 138|) . Thus, it remains to verify (139|) . Since 
for all 771 G N and all h G (0, oo), the Pochhammer like expression T{h + m — l)/(m — 1)! 
appearing on the left hand side in (15^ is bounded below and above by 

1 , f^-^\^b~2 / r(6 + m-l) ,_i /^^-A_b-2 , r^_6-3 



where i^f, := r(6) — 1 — (2^), (l39l) follows by plugging in these lower and upper bounds 
on the left hand side in ( 139|1 and applying afterwards Lemma 16.61 with a replaced by 
a + 6 — 1 > and noting that 



n-l 



^a+6-2 ra°+''"Y ^/logn 



and that 



^-^ [n — rn){n — rn + \)\og^ rn log^n \ \y n 



" 1 ™^a+6-3 /„Q!+b-3 



y ^^ = O 

1^ ^ in — mi In — m -I- 1 i Incr^ m 



_ (n — ?7i) (n — m + 1) log^ m \ log^ n 

by Lemma I6.5[ The proof is complete. D 
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